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We present the dynamic velocity profiles of a Newtonian fluid (glycerol) and a viscoelastic Maxwell 
fluid (CPyCl/NaSal in water) driven by an oscillating pressure gradient in a vertical cylindrical pipe. 
The frequency range explored has been chosen to include the first three resonance peaks of the dy- 
namic permeability of the viscoelastic fluid / pipe system. Three different optical measurement 
techniques have been employed. Laser Doppler Anemometry has been used to measure the mag- 
nitude of the velocity at the centre of the liquid column. Particle Image Velocimetry and Optical 
Deflectometry are used to determine the velocity profiles at the bulk of the liquid column and at 
the liquid-air interface respectively. The velocity measurements in the bulk are in good agreement 
with the theoretical predictions of a linear theory. The results, however, show dramatic differences 
in the dynamic behaviour of Newtonian and viscoelastic fluids, and demonstrate the importance of 
resonance phenomena in viscoelastic fluid flows, biofluids in particular, in confined geometries. 

PACS numbers: 47.50.+d, 47.60.-K, 83.60.Bc 



I. INTRODUCTION 

Coupling between flow and liquid structure 
makes the dynamic response of non-Newtonian 
(complex) fluids much richer than that of Newto- 
nian (simple) fluids [l|,|2j. In particular, depending 
on the relevant time scale of the flow, viscoelastic 
fluids exhibit the dissipative behaviour of ordinary 
viscous liquids and the elastic response of solids. 
Due to their elastic properties, these fluids are po- 
tential candidates to exhibit interesting resonance 
phenomena under different flow conditions. 

In this respect, the response of a viscoelastic fluid 
to an oscillatory pressure gradient has been anal- 
ysed theoretically in some detail. The response, 
measured in terms of the velocity for a given am- 
plitude of the pressure gradient, exceeds that of an 
ordinary fluid by several orders of magnitude at a 
number of resonant frequencies. The remarkable 
enhancement in the dynamic response of the vis- 
coelastic fluid is attributed to a resonant effect due 
to the elastic behaviour of the fluid and the geom- 
etry of the container 0, 0, IU, El • 

This theoretical prediction has been recently con- 
firmed by Laser Doppler Anemometry (LDA) mea- 
surements of velocity at the centre of a fluid column 
driven by an oscillating pressure gradient Q- The 
experiments show that a Newtonian fluid exhibits a 
simple dissipative behaviour, while a Maxwell fluid 
(the simplest viscoelastic fluid) exhibits the reso- 
nant behaviour predicted by the linear theory at 
the expected driving frequencies. While this anal- 



ysis was performed in the central point of the fluid 
column, an exploration of the whole velocity field is 
necessary to ensure that the linear model captures 
the main features of the flow. 

Resonance effects can only be observed when the 
elastic behaviour of the fluid is dominant. This is 
properly characterised by a dimensionless number 
De^> 1, where De (Deborah number) measures the 
relative importance of the relaxation time of the 
fluid to the typical time scale of the flow. More- 
over, at the resonant driving frequencies the system 
response is nearly purely elastic and the dissipative 
effects are very small. 

In the present paper we extend the previous ex- 
perimental investigation in several directions. In 
the first place we present new LDA measurements 
for two Maxwell fluids of different composition, and 
show that the linear theory correctly predicts the 
location of the resonance frequencies in terms of 
material parameters. In addition, measurements at 
different depths along the column centre show the 
attenuation of the velocity as the upper free inter- 
face is approached. The main result is the determi- 
nation of the velocity profiles in the radial direction 
of the cylindrical tube, at several time intervals and 
driving frequencies. We have measured the velocity 
profiles at two different locations: the bulk of the 
liquid column, and the upper liquid-air interface. 
Particle Image Velocimetry (PIV) has been used in 
the former case, and an original technique based 
on Optical Deflectometry (OD) in the latter. The 
PIV results compare satisfactorily with the velocity 



2 



profiles predicted by the linear theory |8| ■ The OD 
results show the attenuation of the velocity due to 
surface tension. 



II. LINEAR THEORY 

In this section we recall the theoretical expression 
of the velocity field in the bulk of a viscoelastic fluid 
contained in a vertical tube^and subjected to an os- 
cillating pressure gradient 0, @ . The tube is sup- 
posed to be infinite in the vertical direction. The 
derivation is based on a linear approximation of the 
hydrodynamic equations and a linearized constitu- 
tive equation for a Maxwell fluid [3|. LiJ. La| • 

Following Ref. 7], the velocity field in Fourier's 
frequency domain reads: 
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In this expression r, z are the radial and vertical 
cylindrical coordinates in a reference system cen- 
tered with the tube axis, u> — 2t:v, where v is 
the driving frequency, t m is the relaxation time of 
the Maxwell fluid, r\ its dynamic viscosity, /3 = 

p [(t m uj) 2 + iwt m ] /(rjtm), a is the cylinder ra- 
dius, Jo is the cylindrical Bessel function of zeroth 
order, and P is the applied pressure in Fourier's 
frequency domain. 

In our experiment, oscillations of the pressure 
gradient are induced by the harmonic motion of a 
piston at the bottom end of the fluid column, with 
amplitude zo (Figure^}. Thus, 
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The velocity profile in the radial direction of the 
tube, at the driving frequency u>, is then given by 
the real part of the following expression: 
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For a Newtonian fluid (t m = 0) there is only one 
node of the velocity profile, at r = a, which ac- 
counts for the non-slip condition at the wall. Re- 
markably, if t m 7^ (Maxwellian fluid) the velocity 
profile may present several nodes. These nodes de- 
fine quiescent points of the flow. For given mate- 
rial and geometrical parameters, the location of the 
nodes depends on the driving frequency, as shown 
in Fig. □ 

As mentioned in the Introduction, the elastic 
properties of the fluid dominate the dynamic re- 
sponse of the system at the resonant frequencies. 
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FIG. 1: Schematic view of the experimental device, in- 
cluding the setups for PIV and defiectometry measure- 
ments. 
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FIG. 2: Diagram showing the location of the quiescent 
flow points along the radial coordinate of the cylinder, 
r, as a function of driving frequency, v, for a Maxwell 
fluid. The dashed horizontal lines separate the fre- 
quency axis in intervals of constant number of quies- 
cent flow points. The parameters p, t m , and r\ used to 
compute the diagram are the ones listed in the text for 
the 60/100 CPyCl/NaSal solution, together with the 
cylinder radius a — 25 mm. 
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In Eq. © this is manifest in that, at the resonant 
frequencies, the oscillation in the centre of the tube 
(maximum velocity) is in phase with the driving. 



III. EXPERIMENTAL DEVICE 

The experimental device is shown in Fig.^ The 
cylindrical container (inner radius a — 25 mm, 
length 500 mm) is made of transparent acrylic. 
In order to avoid optical aberrations this cylin- 
der is placed inside a second recipient of square 
section, made also of transparent acrylic, which is 
filled with glycerol to match the refractive index 
of the acrylic walls. A Teflon piston driven by a 
motor of variable frequency produces harmonic os- 
cillations of the pressure gradient in the liquid col- 
umn. Although the motor controller can operate 
in the frequency range from 1.5 to 200 Hz, our 
experiments have been performed in the reduced 
frequency range from 1.5 to 15 Hz. The piston os- 
cillation amplitude is set to 0.80 ±0.05 mm to keep 
the Reynolds number low, within the accessible fre- 
quency range, for the two fluids studied. Additional 
details of this part of the setup are described in Ref. 


The viscoelastic fluid used in these experiments 
is an aqueous solution of cetylpyridinium chloride 
and sodium salicylate (CPyCl/NaSal). This sur- 
factant solution is known to exhibit the rheologi- 
cal behaviour of a linear Maxwell fluid in a range 
of concentrations |3, Il0| . In most experiments we 
have used a 60/100 concentration, for which the dy- 
namic viscosity 77 = 60 Pa-s, the density p = 1050 
kg/m 3 , and the relaxation time t m = 1.9 s |ll|. In 
DLA measurements we have also used a 40/40 con- 
centration, for which 77 = 30 Pa-s, p — 1005 kg/m 3 , 
and relaxation time around t m = 1.25 s. Our rhe- 
ological characterisation of this last fluid, however, 
suggests that it cannot be fully described by a sin- 
gle relaxation time. As Newtonian fluid we have 
used commercial glycerol, with nominal dynamic 
viscosity 77 = 1 Pa-s and density p — 1250 kg/m 3 . 
These values have been determined at the working 
temperature of (25 ± 0.5)°C. 

All sets of measurements performed on the same 
fluid have been carried out in order of increasing 
frequency. The cylinder has been emptied and re- 
filled with fresh fluid every time that a series of 
measurements at three different driving frequencies 
was complete, i.e. approximately every day. Mea- 
surements in the bulk have been carried out at two 
different heights, corresponding to distances of 6 
and 10 cm from the upper free interface. 

The LDA technique used in the measurements 



presented in the next Section has already been de- 
scribed in detail in Ref. |7j . We present now a brief 
description of the other two techniques (PIV and 
deflectometry). 



A. Particle Image Velocimetry 

The PIV technique provides instant measures of 
the velocity maps in a plane of the flow To 
this purpose, the fluid is seeded with small parti- 
cles. The velocity maps are obtained by a measure 
of the statistical correlation of the displacement of 
the seeding particles in the fluid in a known time 
interval, in this case the time between two consec- 
utive laser pulses. 

Our PIV system contains a two pulsed Nd-YAG 
lasers unit, that includes an optical array to pro- 
duce a laser light sheet in a vertical plane of the 
acrylic cylinder (Fig. A high resolution cam- 
era (Kodak E1.0), perpendicular to the laser light 
sheet, is used to record the digital images. The 
camera records two consecutive frames, one corre- 
sponding to each laser light pulse. The acquisi- 
tion rate is limited by the camera to three pairs 
of images every two seconds (1.5 Hz). A Dan- 
tec FlowMap 1100 processor takes care of the syn- 
chronization between the laser pulses and the cam- 
era trigger. Post-processing of the data, to deter- 
mine velocity maps, is carried out by the Dantec 
FlowMap v5.1 software. Dantec 20-fim polyamid 
spheres were used as seeding particles in the present 
experiments. These particles are small enough to 
follow the flow with minimal drag, but sufficiently 
large to scatter enough light to obtain good particle 
images. 



B. Optical Deflectometry 

Deflectometry is a technique used to study rela- 
tively small deformations in transparent fluid sur- 
faces |l3j . In our case the liquid column serves as 
a variable-thickness lens. A thin transparent plas- 
tic sheet with a regular array (grid spacing d) of 
black dots (0.5 mm diameter) is placed on top of 
the Teflon piston, at the bottom of the acrylic cylin- 
der (Fig.nj. The grid is imaged through the liquid 
column using a CCD camera interfaced to a com- 
puter. Light passing through the liquid column is 
refracted at the liquid-air interface, distorting the 
image of the grid as the interface deforms. Figure 
13 presents an example of the recorded images. 

The displacement of an imaged grid point, 5, is 
related to the local slope of the interface, <j>, in the 
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FIG. 5: Theoretical deformation profiles of the Maxwell 
fluid driven at 2 Hz (dashed line) and at 6.5 Hz (solid 
line) . 
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FIG. 4: Sketch of the path followed by a light ray in 
the system. 



following way (Fig. 0J. According to Snells' refrac- 
tion law, n\9 — n2(j). On the other hand, for a 
column height h the figure shows that tan(0 — 9) = 
S/h. Therefore, the relation between </> and S reads: 



smooth profiles, and the condition of small defor- 
mations is satisfied in a wide range of frequencies. 

Although the deformations are small in this fre- 
quency range, the whole deformation profile de- 
pends strongly on the frequency of driving. This 
can be seen on the theoretical deformation pro- 
files of the viscoelastic fluid at 2 and 6.5 Hz, dis- 
played in Fig. [SJ The two profiles are very differ- 
ent. While the deformation profile at 2 Hz is mono- 
tonic and very smooth, the profile at 6.5 Hz shows a 
more complex behaviour: instead of a single central 
node, it presents two nodes which, in the velocity 
profile, correspond to a maximum at the tube cen- 
tre and a minimum near the tube wall. The max- 
imum deformation between these two nodes corre- 
sponds to an inflection point of the velocity profile, 
very close to the quiescent flow point. 

In our setup the height of the liquid column is a 
function of the radial coordinate r, and periodic in 
time. This height, h(r,t), is simply related to the 
local slope by 4>(r,t) — dh(r,t)/dr. Hence: 



tf> = — — arctan ( ^ ) . (4) 

l~n 2 /ni \hj 

In our setup the height of the liquid column 
(H = 240 mm) is much larger than the local height 
variations produced by the deformation of the in- 
terface (of the order of mm) , and we can safely take 
h ~ H in Eq. ®. 

OD is usually used in thin films. The deflection 
of the points of the grid increases with the depth 
of the layer, making this technique useless for the 
measurement of highly deformed interfaces of thick 
layers. However, the small amplitude of the per- 
turbation in our experiments generates relatively 



h(r,t) = h{r ,t)+ / </>(r,t)dr. (5) 

It is useful to choose for ro a position where the 
interface is motionless, because the reference height 
h(ro,t) is then a constant that can be taken equal 
to 0. We have chosen ro as the radial coordinate 
of the quiescent flow point (in the bulk) closest to 
the cylinder axis. Implicitly, we are assuming that 
the stagnation points predicted by Eq. © for the 
bulk do not change their position as the interface is 
approached. This assumption is confirmed by our 
PIV and OD results (see Sea llVfl . 
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Finally, the velocity profile of the interface can 
be obtained from the time derivative of the local 
height, v(r,t) = dh(r,t)/dt. 

In our experiments the grid spacings used are 
d = 5.0 and 5.5 mm, and the refractive indices 
are n = 1.33407(01) for the CPyCl/NaSal solution 
(measured by Abbe refractometry) , and n = 1.473 
for glycerol (nominal). 



IV. RESULTS AND DISCUSSION 



A. rms velocity at the cylinder axis 
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The amplitude of the velocity field given by Eq. 
I|T|) presents resonance peaks at several resonance 
frequencies. As mentioned in the Introduction, this 
phenomenon was demonstrated experimentally and 
compared to the purely dissipative behaviour of a 
Newtonian fluid in Ref. 7]. These results showed 
that the linear theory gives a good prediction of the 
resonance frequencies but overestimates the ampli- 
tude of the resonance peaks. 

We have repeated this same kind of measure- 
ments in the present work, to check the depen- 
dence of the resonance frequencies on the rheolog- 
ical properties of the Maxwell fluid. To this pur- 
pose the driving frequency to has been rescaled by 
a characteristic time r defined as 



(6) 



where a is the inverse of our Deborah number 
De — t m T)/a p. As it was shown in Ref. |5j, the 
location of the resonance peaks becomes universal 
(independent of fluid parameters and system di- 
mensions) if uj is made dimensionless in the form 

LOT. 

Then, we rewrite Eq. in the form 



(7) 



and plot the dimensionless response function 
£(r, <j)/£(r, 0) as a function of the dimensionless 
driving fequency lot. The results at the centre of 
the tube (r = 0), for the two different concentra- 
tions of the CPyCl/NaSal solution used, are shown 
in Fig. Comparison with the linear theory is 
given by the vertical dashed lines, which give the 
resonance frequencies corresponding to maxima of 
£(r = 0, cj). The Figure shows that the rescaling 
of the frequency axis suggested by the linear the- 
ory leads to a satisfactory reproducibility of the 
resonance frequencies independently of viscoelastic 
fluid parameters. 
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FIG. 6: Dimensionless response function at the centre 
of the cylinder, as a function of dimensionless driving 
frequency, for two concentrations of the CPyCl/NaSal 
solution. The vertical dashed lines give the location of 
the resonance frequencies predicted by a linear theory. 



The magnitude of the response function at the 
resonance peaks for the surfactant solution of con- 
centration 60/100 is considerably lower in the 
present measurements than in the previous ones of 
Ref. [3 • The reason is that the present measure- 
ments have been carried out at 6 cm of the free air- 
liquid interface, while the previous ones had been 
taken at 10 cm. The different measured velocities 
provide a first evidence of the damping influence of 
the free interface on the flow. 

The results presented in the next two sections 
provide measurements of the whole velocity profile, 
instead of measurements at a single point in the 
flow 0. The velocity profiles are obtained at two 
different heights, in the bulk of the fluid, and at the 
fluid-air interface. 



B. Velocity profiles in the bulk 

The following set of figures presents the velocity 
profiles in the bulk of the fluid column, determined 
by PIV measurements, together with the theoret- 
ical profiles given by Eq. © at coincident time- 
phases, for comparison. The profiles have been de- 
termined at the driving frequencies of 2 Hz (Figs.0 
and©, 6.5 Hz (Figs.0 and [TOJ, and 10 Hz (Figs. El 
and El- These values coincide with the first three 
resonance frequencies of our system for the 60/100 
viscoelastic solution, as given by Eq. The time- 
phases have been selected by the criterion that the 
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FIG. 7: Glycerol: PIV results at 2 Hz (dots), measured 
at 6 cm from the upper free interface, and the corre- 
sponding theoretical prediction (solid line) at t — 0.125 
s. 
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FIG. 8: 60/100 CPyCl/NaSal solution: PIV results at 
2 Hz (dots), measured at 6 cm from the upper free 
interface, and the corresponding theoretical prediction 
(solid line) at t = 0.375 s. 



velocity at the tube axis is a maximum. 

The first observation to make is that the in- 
stantaneous velocity profiles of the Maxwell fluid, 
driven at 2 Hz, present a single defined sign of ve- 
locity (single direction of motion) along the whole 
radius of the tube. As the driving frequency is in- 
creased to 6.5 and 10 Hz, however, the instanta- 
neous profiles display a progressively more com- 
plex structure, revealing the presence of annu- 
lar regions within the tube with alternating up- 
ward/downward motion. Notice that this com- 
plexity is inherent to the viscoclastic properties of 
the Maxwell fluid. For the Newtonian fluid (glyc- 
erol) the instantaneous flow in the tube goes all in 
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FIG. 9: Glycerol: PIV results at 6.5 Hz (dots), mea- 
sured at 6 cm from the upper free interface, and the 
corresponding theoretical prediction (solid line) at t = 
0.115 s. 
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FIG. 10: 60/100 CPyCl/NaSal solution: PIV results at 
6.5 Hz (dots), measured at 6 cm from the upper free 
interface, and the corresponding theoretical prediction 
(solid line) at t = 0.038 s. 



the same direction for the three driving frequencies 
tested. 

For the viscoelastic fluid the frontiers between 
consecutive annular regions with alternating signs 
of the velocity do not move. They correspond to 
the stagnation regions of the flow which were al- 
ready discussed in the context of the linear theory. 
The PIV results show that the number of quiescent 
flow points along the radial direction of the tube 
increases with the driving frequency, in agreement 
with the theoretical prediction (Fig. . 

Some of the measured profiles of the viscoelastic 
fluid show regions near the walls with vanishingly 
small velocities and near zero velocity gradients, 
most noticeably the one at 2 Hz shown in Fig. |SJ 
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FIG. 11: Glycerol: PIV results at 10 Hz (dots), mea- 
sured at 6 cm from the upper free interface, and the 
corresponding theoretical prediction (solid line) at t = 
0.075 s. 



FIG. 13: 60/100 CPyCl/NaSal solution: PIV results at 
2 Hz, measured at 6 cm from the upper free interface 
and at different time phases within a driving period. 
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FIG. 12: 60/100 CPyCl/NaSal solution: PIV results at 
10 Hz (dots), measured at 6 cm from the upper free 
interface, and the corresponding theoretical prediction 
(solid line) at t = 0.025 s. 



These profiles are reminiscent of velocity profiles 
obtained for systems that display shear-banding. 
Indeed, the CPyCl/NaSal solution is commonly 
known to show shear-banding but we have 

enough evidence to discard this effect in our experi- 
ments. We did not observe an increase in turbidity 
in the region of the fluid close to the walls, nor 
changes in the local intensity of the scattered light 
that could be attributed to inhomogeneities. Fur- 
thermore, we monitored possible changes on the po- 
larization state of the scattered light, an indication 
of banding . No changes in the optical proper- 
ties were observed for the amplitude and frequency 
range used here. In addition, Fig. demonstrates 
that the velocity near the walls takes values dis- 



tinctly different from zero at different phases within 
an oscillation period. 

Performing PIV measurements at two different 
heights of the liquid column allows studying the in- 
fluence of the upper free interface on the velocity 
profiles. Thus, it is interesting to compare the re- 
sults presented above, which have been performed 
at 6 cm from the upper interface, to the results 
shown in Figs. El an( i which have been per- 
formed at 10 cm from the upper interface. The 
first conclusion to draw from corresponding mea- 
surements at different heights is that the location 
of the quiescent flow points is not affected by the 
presence of the upper free interface. The second 
conclusion is that the magnitude of the velocity 
profile is smaller when the measurement is carried 
out closer to the free interface. The damping ef- 
fect of the free interface, disregarded in the theory, 
originates from the air-liquid surface tension. This 
observation applies equally to the two fluids inves- 
tigated (Newtonian and Maxwellian) . 

The experimental velocity profiles are in reason- 
able agreement with the theoretical ones. The num- 
ber of quiescent flow points is the same, and their 
location very similar. As a general trend, however, 
the theory overestimates the measured velocity. 

A first explanation would be that the theory dis- 
regards nonlinearities, and these tend to limit the 
largest velocity values. Nonlinearities could arise 
from either the hydrodynamic equations or the con- 
stitutive relation of the fluid. In our case, how- 
ever, since the small amplitude of the piston oscilla- 
tions ensures that Re< 10 -4 , the linearized momen- 
tum equation is a very good approximation. On 
the other hand, taking into account the cylindrical 
symmetry of the problem and assuming that the 
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I10I and El reveals that the shear rate actually ex- 
perienced by the fluid in some phase intervals of the 
oscillation is larger Q . In these conditions the vis- 
cosity of the fluid decreases with shear. The theory 
predicts that the dynamic response of the system 
at the resonant frequencies becomes smaller as the 
viscosity is reduced , and thus would support the 
view that the measured velocity profiles are system- 
atically smaller than the theoretical ones for the 
viscoelastic fluid (and not for the newtonian fluid) 
because of shear-thinning. 



FIG. 14: 60/100 CPyCl/NaSal solution: PIV results at 
6.5 Hz (dots), obtained at 10 cm from the upper free 
interface, and the corresponding theoretical prediction 
(solid line) at t = 0.057 s. 
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FIG. 15: 60/100 CPyCl/NaSal solution: PIV results at 
10 Hz (dots), obtained at 10 cm from the upper free 
interface, and the corresponding theoretical prediction 
(solid line) at t = 0.022 s. 



velocity depends only on the radial coordinate (as 
our results confirm to a very good approximation), 
it turns out that the first nonlinear correction to 
the constitutive equation of the fluid cancels out 
exactly. 

The disagreement between theory and experi- 
ment is probably due to shear-thinning of the vis- 
coelastic fluid. As shown in Ref. our fluid is 
properly described by a Maxwell model up to shear 
rates 7 ~ 0.6 s _1 , and experiences shear-thinning 
beyond that value. A close inspection of Figs. |H1 



C. Deflection of the air— liquid interface 

The results obtained by Optical Deflectometry 
are presented on top of Figs. 1161 and!17lf2 Hz), and 
Figs. UHlandUll (6.5 Hz). There are no results at 10 
Hz because, at this driving frequency, the deforma- 
tions of the free interface of the Maxwellian fluid 
are so large that the technique is not applicable. 
Each figure displays measurements at two different 
time-phases. 

The experimental deformation profiles of the in- 
terface at 2 and 6.5 Hz for the viscoelastic fluid can 
be compared to the theoretical ones (computed for 
the bulk) displayed in Fig.JSJ We observe that their 
shapes are fully coincident for each of the two driv- 
ing frequencies. This shows that the quiescent flow 
points of the velocity profiles do not change their 
position along the vertical direction, and thus their 
location in the bulk can be used as reference to 
compute the velocity profiles at the interface from 
the deformation profiles. 

The velocity profiles are presented on the bottom 
of Figs. HI and C3 (2 Hz), and Figs. [H and [H (6.5 
Hz) . The magnitude of the interfacial velocities is 
in all cases much lower than the one measured in 
the bulk of the fluid. This is mostly due to the sta- 
bilizing role of surface tension, as discussed above. 

It is interesting to note that the velocity profiles 
of the interface depend also on the direction of mo- 
tion of the driving piston. The magnitude of the ve- 
locity is systematically lower for positive displace- 
ments of the piston (liquid displacing air) than for 
negative ones (air displacing liquid). We attribute 
this asymmetry to the fact that the large viscosity 
contrast across the interface stabilizes the displace- 
ment of a nearly inviscid fluid (air) by a viscous 
liquid, and destabilizes the opposite displacement. 
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FIG. 16: Glycerol: deflectometry results at 2 Hz, mea- 
sured at two different time-phases within an oscillation. 



FIG. 17: 60/100 CPyCl/NaSal solution: deflectometry 
results at 2 Hz, measured at two different time-phases 
within an oscillation. 



V. CONCLUSIONS 

We have used three different optical techniques 
to characterize the oscillating flow of a viscoelas- 
tic fluid (a solution of CPyCl/NaSal in water) con- 
tained in a vertical cylinder and subjected to an 
oscillating pressure gradient. A Newtonian fluid 
(glycerol) has also been studied for the sake of com- 
parison. 

LDA measurements of the fluid velocity at the 
symmetry axis of the cylinder as a function of driv- 
ing frequency, for two different compositions of the 



surfactant solution, have enabled us to show that 
the frequencies of the resonance peaks can be pre- 
dicted accurately in terms of the fluid rheological 
properties by a simple linear theory. This theory 
neglects inertial effects, and makes use of a linear 
Maxwell model as constitutive relation for the vis- 
coelastic fluid. 

Systematic PIV measurements of the radial ve- 
locity field in the bulk of the fluid have been per- 
formed at three different driving frequencies, close 
to the first three resonance frequencies of the vis- 
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FIG. 18: Glycerol: deflectometry results at 6.5 Hz, 
measured at two different time-phases within an os- 
cillation. 
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FIG. 19: 60/100 CPyCl/NaSal solution: deflectometry 
results at 6.5 Hz, measured at two different time-phases 
within an oscillation. 



coelastic system. While the velocity profile of the 
Newtonian fluid along the radial direction does not 
change sign, this is not the case for the Maxwcllian 
fluid. The profiles measured at 6.5 and 10 Hz 
present regions with alternating signs of the veloc- 
ity, separated by quiescent flow points. The number 
of quiescent flow points increases with the driving 
frequency, revealing the increasing complexity of 
the flow. Measurements within the fluid column 
at two different heights show that these quiescent 
flow points do not shift as one moves along the 



vertical direction, and that their radial location is 
accurately reproduced by the linear theory. 

The presence of a free interface at the top of the 
liquid column has a damping effect on the velocity 
amplitude. This observation is visible, both with 
LDA and PIV, when the results of measurements 
carried out at two different heights within the liquid 
column are compared. 

Optical Deflectometry measurements of the free 
interface confirm that the velocity field is severely 
damped by the surface tension of the air-liquid in- 
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terface, compared to the velocity field within the 
bulk. Interestingly, the deflectometry results show 
also that the oscillations of the velocity field at the 
interface are asymmetric, the profiles correspond- 
ing to positive displacements of the piston having a 
slightly but systematically smaller amplitude than 
those corresponding to negative displacements. We 
attribute this asymmetry to the fact that the up- 
ward motion of the interface (liquid displacing air) 
is stabilized by the viscous pressure gradient in the 
liquid. Recently it has been suggested that the vis- 
cous fingering interfacial instability of a viscoelastic 
fluid in a Hele-Shaw cell will be enhanced by mak- 
ing the flow oscillate at the resonant frequencies 
flq . Our deflectometry results show that the sta- 
bilizing role of the surface tension will reduce this 
effect to some extent. 

These results, as part of the study of resonance 
frequencies in viscoelastics, are potentially relevant 
in biorheology to interpret the frequencies at which 



biofiuids are optimally pumped in living organisms 
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